Setty et al. BMC Systems Biology 201 1, 5:154 
http://www.biomedcentral.eom/1 752-0509/5/1 54 

Systems Biology 



RESEARCH ARTICLE Open Access 



How neurons migrate: a dynamic in-silico model 
of neuronal migration in the developing cortex 

Yaki Setty 1,2 *, Chih-Chun Chen 3 , Maria Secrier 4 , Nikita Skoblov 5 , Dimitrios Kalamatianos 6 and Stephen Emmott 1 




Abstract 

Background: Neuronal migration, the process by which neurons migrate from their place of origin to their final 
position in the brain, is a central process for normal brain development and function. Advances in experimental 
techniques have revealed much about many of the molecular components involved in this process. 
Notwithstanding these advances, how the molecular machinery works together to govern the migration process 
has yet to be fully understood. Here we present a computational model of neuronal migration, in which four key 
molecular entities, Lis 1 , DCX, Reelin and GABA, form a molecular program that mediates the migration process. 

Results: The model simulated the dynamic migration process, consistent with in-vivo observations of morphological, 
cellular and population-level phenomena. Specifically, the model reproduced migration phases, cellular dynamics 
and population distributions that concur with experimental observations in normal neuronal development. We 
tested the model under reduced activity of Lis 1 and DCX and found an aberrant development similar to observations 
in Lis 1 and DCX silencing expression experiments. Analysis of the model gave rise to unforeseen insights that could 
guide future experimental study. Specifically: (1) the model revealed the possibility that under conditions of Lis 1 
reduced expression, neurons experience an oscillatory neuron-glial association prior to the multipolar stage; and (2) 
we hypothesized that observed morphology variations in rats and mice may be explained by a single difference in 
the way that Lis 1 and DCX stimulate bipolar motility. From this we make the following predictions: (1) under reduced 
Lis 1 and enhanced DCX expression, we predict a reduced bipolar migration in rats, and (2) under enhanced DCX 
expression in mice we predict a normal or a higher bipolar migration. 

Conclusions: We present here a system-wide computational model of neuronal migration that integrates theory 
and data within a precise, testable framework. Our model accounts for a range of observable behaviors and affords 
a computational framework to study aspects of neuronal migration as a complex process that is driven by a 
relatively simple molecular program. Analysis of the model generated new hypotheses and yet unobserved 
phenomena that may guide future experimental studies. This paper thus reports a first step toward a 
comprehensive in-silico model of neuronal migration. 



Background 

Neuronal migration is a highly dynamic process that is 
essential for the normal development and function of the 
mammalian brain. The migration process is regulated by 
cell-extrinsic signaling pathways and cell-intrinsic regula- 
tion and implicates numerous molecules that synergisti- 
cally guide neuron motility. The experimental data have 
revealed that the process initiates with proliferation of neu- 
roblasts from the progenitor glial cells at the ventricular 
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zone (VZ). The neuroblasts migrate radially toward the 
subventricular zone (SVZ) along the glial fiber, which 
serves as a their scaffold [1]. On exit from the VZ, neuro- 
blasts adopt a multipolar migration stage in which they dis- 
associate from the glial fiber and migrate independently 
through the intermediate zone (IZ) [1,2]. As the multipolar 
neurons approach the cortical plate (CP), they enter a bipo- 
lar migration stage and re-associate with the glial fibers. On 
entry into the CP zone, the migrating neurons dissociate 
from the fiber and accumulate in a layer above that formed 
by their predecessors, such that layers of progressively 
younger neurons are found from interior to surface [3,4]. 
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Advances in experimental techniques and genetic stu- 
dies have identified key molecular components that are 
implicated in the migration process. Yet, how these com- 
ponents act in concert to mediate neuronal migration at 
the system-level is not fully understood. Key outstanding 
questions that remain unsolved include how neurons 
interpret guidance cues, what are the specific molecular 
mechanisms that guide migration, and how neurons coor- 
dinate interactions between different signaling pathways. 
Theoretical study may afford a platform to explore how 
the different components orchestrate the migration pro- 
cess as a whole. Because of the high degree of molecular 
complexity governing neuronal migration simplification is 
required. This has been constructively demonstrated in 
previous studies that have categorized the molecular enti- 
ties into functional classes [5] or mechanical instructions 
[6]. We have adopted similar simplifications in our model- 
ing study. 

The computational model we present here integrates 
isolated subcomponents of the migration process into a 
dynamic system-level simulation of the process. We 
focused on the interplay between four molecular entities 
that have been found to be crucial for normal migration: 
two extracellular signaling cues, GABA neurotransmitter 
and Reelin protein, and two intracellular regulators, the 
Lissencephalyl (LIS1) and Doublecortin (DCX) proteins. 
We defined the following functionalities for GABA, Reelin, 
Lisl and DCX: (1) GABA neurotransmitter has a chemo- 
tactic role and a chemokinetic effect on migrating neu- 
rons. Experimental observations have found that GABA 
influences migrating cortical neurons [1,7]. We, therefore, 
defined GABA to act as a chemoattractant that provides, 
together with the glial fiber, guidance cues for migrating 
neurons. (2) Reelin is implicated in the maintenance of the 
radial scaffold [8] and neuron polarization as well as in sta- 
bilizing the leading process and promoting the final soma 
translocation [9]. Reelin is secreted by Cajal-Retzius cells 
located at the pial surface and is predominantly expressed 
in the CP [10]. In addition, Reelin was found to induce 
dissociation of neurons from the glial fiber [11]. Therefore, 
in our model Reelin serves as a local stop signal for 
migrating neurons that enter the CP [3,5]. (3) The intra- 
cellular regulators, Lisl and DCX, mediate neuronal polar- 
ization by stabilizing microtubules [12] and are implicated 
in the regulation of leading process branching and 
dynamics [13]. In addition, these regulators play a role in 
the proliferation of neural progenitors [14] and have also 
been attributed with unknown cell-nonautonomous func- 
tions [15]. Silencing expression of Lisl and DCX by RNA 
inhibition has been shown to arrest normal migration and 
to result in aberrant migration with a distinct pattern 
[1,16-18]. Thus, in the model, we specified Lisl and DCX 
as regulators that stimulate and guide neuron motility and 
timing. (4) For modeling purposes we assumed that the 



association with the glial fiber is controlled via the (direct 
or indirect) interaction of the motility regulators, Lisl and 
DCX, with Astrotactin, which was found to mediate neu- 
ron movement along the glial fiber [19]. We are not, how- 
ever, aware of any experimental findings that support or 
refute this assumption. 

Our model incorporates only a subset of the rich genetic 
knowledge. The molecular entities, even in this compact 
subset, have a variety of complex functions and are impli- 
cated in various signaling pathways and mechanisms. 
Nevertheless, we found that the model simulated aspects 
of morphology, distribution and cellular dynamics, consis- 
tent with experimental observations. Furthermore, the 
simulations under reduced Lisl and DCX activity con- 
curred with observations from Lisl and DCX silencing 
experiments using RNA interference. 

The modeling process and the model analysis were 
informative and served as a platform to study the neuronal 
migration at the system-level rather than isolated compo- 
nents. The dynamic nature of the simulation generated 
behaviors that were not explicitly programmed into the 
model. These emergent behaviors are properties of the 
system as a whole and can only be explored as such. The 
simulations give rise to insights that may guide future 
experimental trial. Specifically, the model pointed to the 
possibility that under Lisl reduced expression conditions, 
neurons experience an oscillatory neuron-glial association 
prior to the multipolar stage. Additionally, using the 
model we studied bipolar migration under altered expres- 
sion levels of the key regulators Lisl and DCX in rats and 
mice. The analysis led us to hypothesize that the observed 
morphology variations in these organisms may be 
explained by a single difference in the way the Lisl and 
DCX stimulate bipolar motility. This hypothesis helped us 
to predict the bipolar morphology under yet untested 
expression levels of these regulators. We predict that corti- 
cal migration in rats under reduced expression of Lisl and 
overexpression of DCX will reduce the bipolar migration, 
in contrast with the normal morphology observed in mice. 
Additionally, we suggest that bipolar morphology in mice 
under DCX overexpression will give rise to normal or 
higher bipolar morphology similar to the bipolar morphol- 
ogy observed in rats under similar conditions. 

Results and Discussion 

Model construction and execution 

To construct a model that reproduces the distinct stages 
of the migration process (as illustrated in Figure 1A), we 
specified the interaction scheme shown in Figure IB (see 
Methods for details of model formulation, construction 
and simulation). The model makes a number of simplify- 
ing assumptions: First, interactions occur between five 
extracellular mediators (GABA A , GABA B , ApoER2, 
VLDLR and Astrotactin) that regulate two intracellular 
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Figure 1 Illustration of the neuronal migration process and the simplified set of rules in the in-silico model. A: An illustration of 
neuronal migration (reproduced with permission from [1]). Radial progenitor glial cells (RP; designated by a pink cell with a fiber) proliferate to 
give rise to neuroblasts (PN; designated by a red cell). Neuroblasts migrate guided by the glial fiber at the ventricular zone and subventricular 
zone (VZ/SVZ). Neuroblasts adopt the multipolar stage (MS; designated by a green star-like cell) when entering the intermediate zone (IZ) and 
migrate independently from the glial fiber (on both X and Y directions). In the cortical plate (CP) neurons adopt bipolar state (designated in 
yellow), associate with the glial fiber and accumulate on top of the pile. The histology images correspond to the various migration stages of the 
migrating neurons: neuroblasts (bottom-left), multipolar stage neurons (right) and bipolar neurons (top-left). B: Description of the interaction 
scheme in the model: GABA A detects low concentrations of GABA neurotransmitter in the cortex, but maintains its activity for short periods. 
GABA A in its active form promotes the activity of Lis 1 . GABA B requires higher concentrations of the neurotransmitter. Active GABA B promotes 
both Lis1 and DCX independently. Active Lis 1 and active DCX promote Astrotactin adhesion with the glial fiber. The two Reelin Receptors 
(VLDLR and ApoER2) act as one unit that detects Reelin to reduce the ability of Astrotactin to adhere with the glial fiber one activated. 

V ) 



regulators (Lisl, DCX). A decreasing GABA concentra- 
tion gradient exists from the CP toward the VZ [1,7] 
reflected by a trend in Reelin activity over the CP [3,20]. 
Our second assumption is that once the GABA A receptor 
of a neuroblast detects GABA, it promotes Lisl activity 
to maintain the glial fiber association [21]. This assump- 
tion is based on previously published possible relation- 
ships between GABA and Lisl in C. elegans [22], and 
in-vitro experiments that suggested that GABA induces 
dissociated neurons to migrate from the VZ to the IZ [2] 
via GABA A receptor, which possibly plays role in neuron 
movement [7]. Third, we specify that GABA A activity is 
reduced to a level below that required to promote Lisl. 
Consequently, neurons enter the multipolar stage and 
disassociate from the glial fiber to migrate randomly in 
both X and Y directions [1,23]. Fourth, when GABA B 
detects higher concentrations of GABA, it activates Lisl 
and DCX to re-associate with the fiber and to enter the 
bipolar stage. This is based on observations that GABA B 
receptors promote migration from the IZ to the CP [1,7], 
and on a suggested role for DCX in stabilizing neurons 
association with the glial fiber [2]. Finally, bipolar neu- 
rons migrate radially until the Reelin receptors (i.e., 
ApoER2 and VLDLR) detect and bind to Reelin (i.e., the 
stop signal) and inhibit neuron's adherence ability to the 
glial fiber. In turn, neurons dissociate from the glial fiber 
and accumulate on the next layer in the CP [4,20]. To 
simplify the complex mechanism by which Lisl, DCX 



and Reelin regulate the association with the fiber, we 
defined for modeling purposes that these molecular com- 
ponents mediate the adherence ability of neurons to the 
fiber through interaction with the Astrotactin protein 
that mediates the neuron-glial adhesion ligand (although 
no experimental evidence supports this assumption) [19]. 

The model is seeded by 23 radial progenitor glial cells 
in the simulated cortical cross-section. The progenitor 
glial cells proliferate and give rise to neuroblasts. Each of 
the neurons executes the same molecular program to 
determine its migration behavior over time. For example, 
a neuroblast at the SVZ activates GABA A receptors in 
response to GABA concentrations in its vicinity (see 
detailed description in Methods). 

The model recapitulates complex system-level neuronal 
migration 

Model executions provided a dynamic representation of 
neuronal migration in the cortex (Figure 2A and the 
movie in Additional File 1). Neuroblasts (red) proliferated 
from the glial mother cells (pink) and migrated along the 
glial fiber in the VZ/SVZ radially toward the IZ. On entry 
to the IZ, neurons adopt the multipolar stage and migrate 
randomly both horizontally and vertically independent 
from the glial fiber (green). Finally, neurons adopted a 
bipolar state (yellow) and re-associated to a glial fiber. 
The neurons maintained glial-guided radial migration 
when locating from the IZ to the CP. Having entered the 
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Figure 2 Normal neuronal migration in the mammalian cortex. A. Snapshots of the neuronal migration simulation. Radial progenitor glial 
cells (pink) proliferate asymmetrically giving rise to neuroblasts (red). Neuroblasts migrate radially on the glial fiber. Neuroblasts adopt a 
multipolar stage (green) and migrate independently from the glial fiber. Multipolar neurons re-associate with the glial fiber, enter the bipolar 
migration (yellow) and accumulate on the pile surface. B. Histology of the mouse cortex at embryonic day 16 (reproduced with permission from 
[53]). 



CP, they continued vertical bipolar migration before dis- 
sociating from the fiber in response to Reelin and subse- 
quently accumulating in layers at the CP. 

Simulation results and experimental observations of 
neuronal population distribution were in general agree- 
ment. The model, consistent with histology (Figure 2B), 
generated a narrow band of progenitor glial cells at the 
VZ/SVZ boundary, a sparse occupancy of neurons in the 
IZ and neuron accumulation at the cortical surface (CP). 
The multipolar population dominated the earlier stages 
of migration (blue) and gradually adopted a bipolar fate 
(red) (Figure 3A). Similarly, the occupancy of neurons in 
the different zones of the simulated and experimental 
cortex cross-section showed a gradual migration over 
time. At the early developmental stages the majority of 
neurons positioned at the VZ/SVZ (blue). As the simula- 
tion advances the neurons rapidly move to the IZ (green) 
and then accumulate at the CP (red) (Figure 3B). At the 
end of the simulation 77% of the neuronal population 
was bipolar, 14% multipolar and 9% radial glial cells. The 
cortex occupancy in the simulated cross-section showed 
9% of cells occupying the VZ/SVZ zones, 11% the IZ and 
80% the CP. Experimental observations report a similar 
population distribution in the mouse cortex [24]: 75% 
bipolar neurons, 10% multipolar and 10% glial cell (5% 
unclassified), with zone occupancies of 7%, 15% and 78% 
at the ZV/SVZ, IZ and CP zones respectively. 

Neuron behavior differs between different migration 
stages. Multipolar neurons actively explore their environ- 
ment for directional cues and migrate independently 
from the glial fiber, moving randomly in the X and Y 
directions [2,25]. In contrast, bipolar migration is glial- 
guided and is characterized by a radial movement biased 
in the Y direction. Previous studies [17,25] suggested that 
it is possible to distinguish between multipolar neurons 



and bipolar neurons based on the path they follow and 
their respective velocities. To test the relationship 
between migration stage and path formation, we applied 
a previously employed mathematical technique that dis- 
tinguishes between multipolar and bipolar neurons by 
path formation in-vitro [25]. The analysis revealed a 
mean square error of less than 10% between the classifi- 
cations of neurons in the model using in-vivo and in- 
silico parameters (Figure 3C). Additionally, we analyzed 
relationships between the velocity of the neurons and 
their migration stage by plotting migration stage as a 
function of the horizontal and vertical components of the 
migration velocity (Figure 3D). Neurons in the multipolar 
stage were found to have a low velocity in both directions 
(blue), while neurons in the bipolar stage had typically a 
higher, predominantly Y-directional, velocity (red). This 
analysis concurs with experimental findings that multipo- 
lar neurons remain closer to their origin, whereas bipolar 
neurons are found at a significant distance from their ori- 
gin [17,25]. 

Abnormal GABA-Lisl interaction anticipates aberrant 
migration as observed in Lisl RNAi studies 

Lisl regulator has a complex function in the migration 
and plays a role in various signaling pathways and 
mechanisms. One function that is attributed to Lisl is sta- 
bilizing neuron association with the glial fiber [2]. Normal 
migration in Lisl RNAi mice is arrested in its early stages 
[1,2,21] and displays accumulation of multipolar neurons 
at the SVZ and lower IZ [21]. In our model, Lisl acts as a 
motogen that stimulates and guides the motility of neu- 
rons along the glial fiber in response to GABA concentra- 
tions. Thus, in Lisl reduced activity simulations, 
neuroblasts that detect GABA concentrations did not 
express Lisl normally. Under these conditions, neuroblasts 
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Figure 3 Cellular and population dynamics in neural migration simulations. A. Distribution of multipolar (blue) and bipolar (red) neurons 
over time (averaged over 10 different simulations). B. Zone occupancy as function of time (VZ/SVZ blue), IZ (green) and CP (red) (averaged over 
10 different simulations). C: Migration stage over time as determined by the path formed by neuronal movement (Solid line: mathematical 
analysis of the migration data using in-vivo parameters (from [24]) multipolar population (blue) bipolar population (red). Circles: mathematical 
analysis of the migration data using in-silico parameters: multipolar population (green) and bipolar population (purple). D: Migration stage as a 
function of velocity in X and Y directions, multipolar (blue) and bipolar (red). 



expressed Lisl only at the random basal level embedded in 
the molecular program. 

In simulations under reduced Lisl activity, neuroblasts, 
once proliferated from the mother progenitor glial cell, 
entered the multipolar stage shortly after proliferation and 
dissociated from the glial fiber, rather than progressing nor- 
mally along the fiber. Thus, reduced Lisl activity shortens 
the duration over which neuroblasts remain associated with 
the glial fiber, and hastens multipolar migration. Conse- 
quently, the neuron population accumulated at the SVZ/ 
VZ and the lower IZ and rarely migrated to the CP (Figure 
4A and the movie in Additional File 2). This migration pat- 
tern that emerged in Lisl reduced activity simulations con- 
curs with the morphology observed in Lisl RNAi in mice 
[21] in which neurons accumulate at the lower parts of the 
brain section closer to the VZ/SVZ (Figure 4B). 



The population distribution in the Lisl reduced activity 
simulations showed predominant multipolar cells that 
rarely adopt a bipolar state or enter the CP (Figure 4C, 
D). At the final stage of the simulation, the population 
consisted of 2% bipolar, 89% multipolar and 9% glial 
cells, of which 55% occupy the VZ/SVZ area, 43% the IZ, 
and 2% the CP. These results concur with experimentally 
measured distribution of the neuron population as 4% 
bipolar, 70% multipolar, 10% glial cell, (16% unclassified), 
and zone occupancies of 70%, 20% and 10% in VZ/SVZ, 
IZ and CP, respectively [24]. 

Under simulated conditions of reduced Lisl activity, neu- 
rons lacking normal Lisl-activity responses to GABA 
showed altered duration and timing of association with the 
glial fiber, with subsequent cell accumulation at the SVZ 
and lower IZ areas. The emergent pattern in these 
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Figure 4 Lisl reduced activity experiment. A. In-silico migration in Lis 1 -defective simulations. B. Histology of neuronal migration in the cortex of 
a Lisl RNAi in mice (reproduced with permission from [24]). C. Distribution of multipolar (blue) and bipolar (red) neurons over time (averaged over 
10 different simulations). D. Zone occupancy as function of time (VZ/SVZ blue), IZ (green) and CP (red) (averaged over 10 different simulations). 



simulations points to a possible direct role of Lisl -GAB A 
interaction in the temporal progression of neuroblasts in 
the early stages of neuronal migration. This hypothesis 
could be tested experimentally by increasing GABA con- 
centration in the cortex of animals with suppressed Lisl 
activity. The higher GABA concentrations may boost Lisl 
expression and restore normal association, possibly repair- 
ing the migration pattern in the reduced-expression strains. 

Lisl reduced activity simulations predict an oscillatory 
neuron-glial association prior to the multipolar stage 

Analysis of the neuron-glial association in simulations of 
reduced Lisl activity revealed an intriguing unforeseen 
phenomenon. We found that neuroblasts repeatedly 
attempt to associate with the glial fiber, displaying an 
oscillatory association and dissociation with the fiber, 
especially during the early stages of the neuron lifecycle. 
The oscillatory behavior began immediately following 
neuronal proliferation and lasted until completion of the 
multipolar migration stage. Thus the model, although 
not explicitly programmed to do so, produced a delay in 
the neuronal multipolar migration stage process (Figure 
5 and the movie in Additional File 3). 

The frequency and the duration of the oscillatory delay 
varied between individual neurons: 31% of the neurons 
began multipolar migration immediately after prolifera- 
tion (Figure 5A) 58% had a delay of one embryonic day 
or less (Figure 5B) and 11% maintained this migratory 



behavior for longer than one embryonic day (Figure 5C), 
with 4% exhibiting this behavior over the entire simulated 
period (Figure 5D). Neither the oscillatory frequency nor 
the duration, were explicitly defined in the model. Rather, 
they emerged from the specified role of Lisl and are pro- 
moted by the basal activity of the regulator. 

The periodic association-disassociation phenomenon 
and the delay it forced in the multipolar migration 
could arise from contradicting instructions under 
reduced Lisl activity. The neuroblasts received a weak 
signal to 'associate' from the basal level of Lisl activity, 
while being instructed to 'migrate randomly' by the 
modified Lisl activity. 

We propose that the abnormal motility and defective 
neuronal division observed in Lisl RNAi [21] emerge 
from a delay in the multipolar stage. This finding sup- 
ports a hypothesis by which Lisl plays a significant role 
in neuroblast development and cell motility [26]. These 
results point to a future experimental study to test the 
effect of the combination of timing and consequent spa- 
tial position in individual neuroblasts under reduced 
Lisl activity conditions. 

Simulations under reduced DCX activity anticipate the 
aberrant migration observed in DCX RNAi rats and not 
the normal migration observed in DCX RNAi mice 

DCX function is highly complex and implicates various 
mechanisms. One of its possible functions is that it acts 
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Figure 5 Delay in multipolar migration in Lis 1 -defective neurons. Oscillatory association-dissociation behavior of neurons from the glial fiber 
in four representative individual neurons. A. A cell that immediately adopts multipolar migration. B. Short delay. C. Long delay. D. Oscillatory 
association over the entire simulated period. 



as a motility regulator that increases glial guided migra- 
tion [2,27-30]. DCX RNAi rats show disruptions in 
radial migration, forming a subcortical band consisting 
of multipolar neurons in the IZ [31,32]. In our model, 
DCX acts as a motility regulator that stimulates re-asso- 
ciation with the glial fiber in multipolar neurons. When 
neuronal ability to activate DCX is reduced to the basal 
level, re-association of multipolar neurons with the glial 
fiber is less likely to occur. Consequently, in simulations, 
the bipolar migration stage was largely halted, neurons 
did not advance from the IZ to the CP and reduced 
accumulation in the CP occured. Simulations under 
reduced DCX activity displayed accumulation of multi- 
polar neurons at the upper end of the IZ closer to the 
CP. A minority of the neurons reacted to the random 
basal activity level and adopted the bipolar state, allow- 
ing them to enter the CP (Figure 6A and the movie in 
Additional File 4). The emergent pattern in reduced 
DCX activity simulations is consistent with observations 



that DCX RNAi rats develop a subcortical band of mul- 
tipolar neurons [1,31,32] (Figure 6B). 

Over the simulated period, multipolar neurons concen- 
trated predominantly in the IZ, with only a small propor- 
tion (7%) adopting the bipolar state and entering the CP 
(Figure 6C-D). By the end of the simulation, the neuron 
occupancy in the different zones was in agreement with 
measurements from DCX RNAi rats [27] (numerical 
occupancies of 9%, 55% and 36% vs. experimentally 
observed occupancies of 10%, 60%, 30%, in the VZ/SVZ, 
IZ, and CP zones respectively). The final simulated neu- 
ron population comprised 9% glial cells, 82% multipolar 
stage cells and 9% bipolar stage cells. Though direct 
experimental measurements are not available for com- 
parison, multipolar neuron population was increased and 
bipolar population was decreased in DXC RNAi rats, 
[30,33,34] supporting the simulation's results. 

These results might explain the occurrence of DCX 
re-expression in DCX RNAi rats, with the resumption 
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Figure 6 DCX reduced activity experiment. A. In-silico migration in DCX-defective simulations. B. Histology of neuronal migration in the cortex 
of rats in DCX RNAi (reproduced with permission from [31]). C. Distribution of multipolar (blue) and bipolar (red) neurons over time (averaged over 
10 different simulations). D. Zone occupancy as function of time (VZ/SVZ blue), IZ (green) and CP (red) (averaged over 10 different simulations). 
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of normal bipolar migration [35]. When DCX is re- 
expressed, neurons are able to re-associate to the glial 
fiber, resulting in a normal bipolar morphology. Simi- 
larly, a highly bipolar morphology in DCX over-expres- 
sion rats [27,31] could be explained by an intensive re- 
association of neurons that elevates the bipolar state. 

A single change in the way Lisl and DCX guide motility 
may explain phenotypic variations between rats and mice 

Studies investigating neuronal migration in Lisl and 
DCX RNAi have typically focused on two model organ- 
isms, the rat and the mouse. Lisl RNAi induced abnor- 
mal neuronal migration in both organisms [1]. In 
contrast, neuronal migration in DCX-RNAi animals 
showed a distinct phenotypic pattern: DCX RNAi in rats 
disrupted radial migration in the cortex [27,31], whereas 
DCX RNAi mice developed normally [31]. One possible 
hypothesis is that at the bipolar stage the mutual activity 
of Lisl and DCX is different in the rat and the mouse. 
This hypothesis is supported by observations in Lisl 
RNAi mice in which migration was largely restored 
under DCX over-expression [35], and concurs with pre- 
viously suggested reciprocal functions and possible 
interaction between Lisl and DCX [16,32,36]. 

In our model, reduced DCX activity generated an 
abnormal bipolar morphology consistent with that 
observed in DCX RNAi rats, rather than the normal 



morphology observed in the DCX RNAi mice [31]. DCX 
activity in the model is independent of Lisl since both 
regulators are triggered autonomously in multipolar neu- 
rons. Thus Lisl functioned normally in DCX reduced 
activity conditions, suggesting that Lisl activity alone is 
not sufficient to promote normal bipolar migration. 

These experimental observations combined with the 
model behavior led us to hypothesize that the rat cortex 
requires expression of both DCX and Lisl for neurons 
to enter the bipolar migration stage. In contrast, bipolar 
morphology in the mouse cortex requires that neurons 
either normally express Lisl or over-express DCX 
(Table 1). From this we predict that since the rat cortex 
requires expression of both Lisl and DCX, higher 
expression of DCX in Lisl RNAi rats will not restore 
normal migration as observed for Lisl -defective mice 
[36]. In addition, since Lisl activity is sufficient to pro- 
mote bipolar migration in mice, over-expression of DCX 
in this organism will have no effect on bipolar migra- 
tion. This prediction allows for the possibility that DCX 
over-expression in mice enhances bipolar morphology 
as has been observed in rats [30,33,34] (Table 1). 

Conclusions 

Our computational model of neuronal migration during 
brain development integrates a set of molecular sub- 
components that have been studied in isolation. In the 
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Table 1 Bipolar migration under varying expression levels of Lisl and DCX regulators 



Lisl DCX Bipolar Migration (mouse cortex) Bipolar Migration (rat cortex) 



Reduced 


Reduced 


Reduced 


Abnormal 


Reduced 


Normal 


Reduced 


Abnormal 


Reduced 


Over 


Normal 


Abnormal (*) 


Normal 


Reduced 


Normal 


Abnormal 


Normal 


Normal 


Normal 


Normal 


Normal 


Over 


Normal/Enhanced (*) 


Enhanced 



Bipolar morphology under different expression levels of DCX and Lisl in the mouse and rat cortex. The asterisk symbol indicates predicted morphologies. 



model, cells are represented as executable objects whose 
migration is determined by explicit interactions between 
molecular entities defined and specified as a molecular 
'program'. Overall simulated migration is thus an out- 
come of the collective molecular-level interactions in 
each object In-silico experiments were conducted by 
changing parameters, enabling hypotheses and predic- 
tions about the process to be tested. This model is a 
first step in utilizing in-silico models and computational 
approaches for neuronal migration research. 

The model described here focuses on the interplay 
between two major intracellular regulators, Lisl and 
DCX, the glial fiber, and two extracellular guidance and 
motility cues, GABA and Reelin. While these are a subset 
of key molecular components, they nevertheless repre- 
sent functional classes of mechanisms implicated in neu- 
ronal migration: guidance cues (GABA and the glial 
fiber), motility regulators (Lisl and DCX), and stop sig- 
nals (Reelin). 

The model generated complex dynamics and behaviors 
of the migration process that are in qualitative agreement 
with the migratory process of cortical neurons in normal 
and abnormal conditions. Specifically, our model repro- 
duced the distinct migration stages in the known cortex 
zones [37], bipolar and multipolar cell states over time 
[1,2], cell population dynamics [17,38], and specific disrup- 
tions of the migration induced by manipulations of the key 
regulators, Lisl and DCX [24,27]. Furthermore, simulated 
condition of reduced Lisl activity revealed that neuroblasts 
repeatedly associate and dissociate with the glial fiber 
before initiating a glial-independent multipolar migration. 
In addition, we hypothesized that a single change in the 
way neurons interpret the guidance cues can explain phe- 
notypic variation between different organisms. These 
hypotheses could guide future experimental study. 

We simplified the model for this study, deliberately 
omitting several aspects of migration. Obvious next steps 
would be to undertake more comprehensive modeling of 
the role of Lisl, DCX, GABA and Reelin, for example by 
a gradual activity. Similarly, the model could be extended 
to include (i) cell density within the different regions, (ii) 
neuron-neuron interaction, (iii) neuron-glial interaction, 
(iv) retarding forces encountered by migrating neurons, 



and (v) geometry of migrating neurons. Further elabora- 
tions might incorporate the role of migration factors we 
have excluded from the current study; for example, intra- 
cellular effectors known to influence neuronal migration 
such as CDK5 [39], signaling factors (e.g., Semaphorin- 
3A) [40] and calcium signaling [41]. Similarly, for the 
receptors implicated in migration, a possible refinement 
could incorporate in the model the effects of neuropilin- 
1 [40] and connexins [42]. 

More broadly, the neuronal migration model presented 
here highlights the potential for integrating theory, com- 
putational approaches and experimental data in neu- 
roscience and in biology more generally [43-45]. Indeed, 
increasing focus is given to in-silico approaches for sys- 
tem-level analysis of the dynamics of complex biological 
systems (In-silico study of pancreatic morphogenesis and 
differentiation. Submitted; and [46-51]). We envision that 
in the long run, computational models will become impor- 
tant tools for understanding biological systems, for the 
identification of precise causes of abnormal development, 
and possibly disclose ways to restore normal development 
and function. 

Methods 

Setting up a computational platform 

To compile the specifications of the model into an execu- 
table program (the Java source code of the model can be 
found in Additional File 5), we used MASON (Multi- 
Agent Simulator of Neighborhoods, http://cs.gmu.edu/ 
-eclab/projects/mason/) coupled with Dynamic Data Dis- 
play, a visualization tool developed by Microsoft Research 
(http://dynamicdatadisplay.codeplex.com/). MASON is an 
extendable, discrete event multi-agent simulation toolkit 
designed for multi-agent simulation tasks. This enables 
the model to simultaneously simulate operations and 
interactions of multiple agents, each operating a program 
using states and events that cause transitions between 
states, allowing the system to simultaneously drive several 
orthogonal states (designed in statecharts [52], using prin- 
ciples previously applied in [46]). Dynamic Data Display 
provides visualization of both current state and dynamic 
change from the previous state. The application uses linear 
interpolation algorithms to build intermediate values for 
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the discrete states to visualize changes of the system over 
time. We implemented the platform in C# and used WPF 
library for User Interface purposes and data binding 
between the datasets and the animated objects. Visualiza- 
tion was coded using Microsoft DirectX 9 technology 
(SlimDX DirectX 9 managed wrapper; slimdx.org). 

Specifications of the neuronal migration model 

A cortex cross-section is defined in the model as a 2D 
grid that overlays the agents and represents a cross sec- 
tion of the cortex. Glial fibers over the different cortical 
zones are represented by twenty-five vertically spanning 
routes (see illustration in Figure 7). The ventricular/sub- 
ventricular, the intermediate zone and the cortical plate, 
each span a third of the cortical cross-section. The grid 
specifies concentrations for GABA and Reelin in a way 
that matches their concentrations in real cortex; that is, 
Reelin concentration is specified in the CP, and GABA 
concentration is specified as a gradient from the CP (top 
area of the grid) to the SVZ (bottom area of the grid). As 
an approximation, we assumed a linear gradient for 
GABA and Reelin degradation. A neuron is represented 
by an agent whose activity is governed by the specifica- 
tions described diagrammatically in Figure 8. The neuron 
agent consists of three different objects, Cell, Membrane 
and Effectors. The Cell object comprises two compo- 
nents, Migration Stage and Proliferation. The Migration 



> 
1/1 

-< 
INJ 
> 



< 
CD 
< 



iii 




V 7 



Figure 7 Illustration of the simulated cortical cross-section Glial 
cells, whose nucleus positions lie at the bottom of the cross-section 
and their fiber is stretched to the surface of the cortex (of the cross- 
section). Reelin and GABA factors are present at densities that lie in 
a gradient from top to bottom of the cross-section. 



Stage component specifies possible migration stages for 
the agent. Two cell types, glial cells and migrating neu- 
rons, are distinguished by their specific roles in the simu- 
lation. Glial cells remain in the Glia state and proliferate, 
whereas neuron agents change their migration mode and 
evolve through the different developmental stages - Neu- 
roblasts, Multipolar Stage and Bipolar Stage. The Prolif- 
eration component specifies five cell cycle phases GO, Gl, 
S, G2 and Mitosis. Once a glial agent enters the Mitosis 
stage, a new agent is created and enters the Neuroblast 
state. The Effectors object serves as a cell- intrinsic regu- 
lation unit that specifies behaviors of two effectors: Lisl, 
DCX, each of which can be in an Active or Inactive state. 
The Membrane object serves as a cell-extrinsic signaling 
unit for the agent and consists of receptors that sense the 
cell's vicinity in the grid, and a motility unit that controls 
cell motion. The model has three components for the 
GABA A , GABA B and Reelin Receptors (which account 
for both VLDLR and ApoER2). Each receptor specifies 
activity in a binary fashion defined by an Unbound or 
Bound state. Similarly, the membrane specifies two states 
for the Astrotactin protein, Associated and Not Asso- 
ciated (with the glial fiber). The Motility component 
senses the close vicinity of the cells, continuously seeks 
possible moves and places the agent in a new position if 
applicable. The motion differs between agents in the 
states. Glial cells move vertically in four positions 
between the VZ and the SVZ during a cell cycle. Neuro- 
blasts with glial-cell-associated Astrotactin migrate along 
the fiber. Multipolar neurons migrate randomly, whereas 
bipolar agents migrate vertically, provided that Reelin is 
not present. When Reelin is detected, the agent sets 
Astrotactin to Not Associated and calculates the nearest 
position at minimal distance from the IZ. To incorporate 
basal activity of regulators and receptors resulting from 
spontaneous intrinsic expression, we added noise to the 
transitions (i.e., a random element for transitions 
between the states). A cell may move or make a transi- 
tion (e.g., from inactive to active) with a probability of 
0.01. These transitions occur spontaneously in the 
absence of designated triggering events. We verified that 
the random element itself could not stimulate proper 
migration. Lisl- and DCX-defective simulations were 
implemented by erasing the Promote transition in the 
corresponding component; no additional constraints 
were forced upon the model. See the supplementary 
information in Additional File 6 for the result of a control 
experiment and the designs of the reduced activity in- 
silico experiments. 

Dynamic interactions in the simulation at run-time 

Once the simulation executes, the cortex cross-section 
with 23 glial cells is initiated with 23 agents whose 
Migration Stage component is set to Glia and which are 
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Figure 8 Design of a neuron agent in the model. Specifications of the model of a neuron as an agent in the computational model. Each 
agent consists of three objects, Effectors for intracellular regulation, Membrane for extra cellular signaling and Cell for developmental stages and 
proliferation. Each object specifies possible states for the element it covers. At run time, each orthogonal component (separated in dashed line) 
can be in one active state at a time. 



positioned at the bottom of the grid (VZ). The Prolifera- 
tion component of these cells continuously moves over 
the cell cycle stages and concurrently, the motion unit 
moves the agent vertically over the four bottom fiber 
indexes in the grid. When the agent enters the Mitosis 
stage, a new agent, whose state is set to Neuroblast is 
created. As the proliferation terminates, the Motility 
Unit places the glial cell in the initial position for GO 
phase. Initially, the Effectors of the new agent are in an 
Inactive state. GABA A , GABA B and Reelin Receptors are 
set to Unbound and the Astrotactin is in the Associated 
state. As the simulation advances, the Motion Unit of 
the agent directs its migration radially over the glial 
fiber indexes in the grid. As the neuron moves along the 
glial fiber, it detects GABA concentration and sets the 
state of the GABA A receptor to Bound. In turn, the cell 
promotes the Lisl effector, which enters the Active state 
to maintain normal association with the fiber. On entry 
into the intermediate zone, Lisl activity is reduced as 
the receptor loses its activity and moves its state to 
Unbound. In turn, the cell sets its state to Multipolar, 
the Astrotactin state to Not Associated and the motility 
unit directs random movement over the 2D space. 
When the GABA B receptor detects the higher GABA 
concentrations (i.e., the values in the grid beyond a 



predefined threshold) in the middle grid area (IZ zone), 
it sets its state to Bound and consequently sets the state 
of DCX and Lisl to Active. The neuron then enters the 
Bipolar state, its Astrotactin becomes Associated and 
the Motility Unit directs a radial motion. Finally, when 
the Reelin receptors component detects Reelin in the 
cortex cross-section, its state is set to Bound, the Astro- 
tactin moves into the Not Associated state and the 
Motility Unit searches for an optimal location to 
accumulate. 

Mathematical analysis of the simulation outputs 

To determine multipolar and bipolar stage of agents 
from the path formed at a specific time point, we 
employ the motion directionality analysis in [25] as fol- 
lows: 



where d = J{x{t k+1 ) - x{t k )f + (y(t fe+ i) - y{t k )f is the 
net displacement during a time subdivision (t k ), fe+i) — 
(0, t n ) for which the neuron preserves the same migra- 
tion stage, where t n is the final time point of the simula- 
tion, 0 < k < n, and / = T(< Ux> 2 + < Uf> 2 ) (for a time 



Setty et al. BMC Systems Biology 201 1, 5:154 
http://www.biomedcentral.eom/1 752-0509/5/1 54 



Page 12 of 13 



interval T) is the length of the neuron path. The index J 
scores values in the range between 0 and 1 to specify 
directionality of an agent motion (0 for complete ran- 
domness and 1 for movement along a straight line). 
This parameter quantifies the amount of entropy in the 
movement of a neuron by comparing the length of the 
total path of the cell between two different time points 
and the straight distance between them. It assesses 
whether the movement of the cell lean toward determi- 
nistic or random. This parameter can therefore be used 
to distinguish between the more directed bipolar move- 
ment and the more random multipolar migration. 
Therefore, the higher the /-index score, the more likely 
the neuron is bipolar, and vice versa; the lower the 
score, the more likely the neuron is multipolar. The /- 
index threshold that distinguishes multipolar from bipo- 
lar was determined by an optimization algorithm aimed 
at maximizing agreement with the multipolar population 
data at three time points: 1/3, 2/3 and end of simula- 
tion. The optimization approach uses a simulated 
annealing algorithm composed of successive trials, each 
of which sets new thresholds as the arithmetic mean of 
the two optimal previously computed thresholds, with 
noise insertion to avoid local maxima. The optimization 
yielded a range of threshold values 0.2 to 0.4 inclusive, 
from which we set the threshold to be 0.25 as in [25]. 
The velocities of neurons are computed for each time 

dx 

interval (t k ), (t k+ i) <= (0, t n ), as — i.e., the ratio of the 

dt 

path length for the defined time interval to the length of 
the respective time interval: 



n x {t k+ i) = 



" x t k 



tk+i - h 



and n y (t k+ i) = 



Yt k+ i Yt k 
tk+i - h 



With n x and n y being the horizontal and vertical velo- 
city components and %t i} Yu being the positions of the 
neuron at time t b 0 < i < n ; t n is the final time point of 
the simulation and 0 < k < n. 

Additional material 



Additional file 1: Recorded simulation of neuronal migration under 
normal conditions. 

Additional file 2: Recorded simulation of neuronal migration under 
reduced Lisl activity. 

Additional file 3: Oscillatory association-dissociation behavior of 
four representative neurons. Neurons are labeled green when they are 
disassociated from the glial fiber and change their color to red upon 
association. 

Additional file 4: Recorded simulation of neuronal migration under 
reduced DCX activity. 

Additional file 5: A compressed file of the source code of the 
model. Written in Java and compiled using Mason (see Material 
and Methods). 

Additional file 6: Supplementary information: a control experiments 
and the design of the reduced activity in-silico experiments. 
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